Entanglement of thermal scalar fields in a compact space 
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Using the thermal Green's function approach we propose a general method to investigate entan- 
glement of the vacuum state or thermal ground states in an arbitrary dimensional space-time. As 
' an application we show quantum separability of the massive thermal scalar field vacuum in the 1+1 

dimensional cylindrical space-time. Separability is demonstrated using the positive partial transpose 
criterion for effective two-mode Gaussian states of collective operators. In this case, for all mass 

■ and temperature values entanglement is absent between the collective operators. 
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PACS numbers: 03.67.Mn, 03.67.-a, 03.65.Ud, 71.10.Ca 

' Entanglement is now treated as a physical quantity, as well as valuable resource allowing various quantum infor- 
mation processing. This brings a renewal of interest in studying the non-locality and entanglement in many particle 
' systems such as Bose Einstein condensations [l[ , fermion systems [U, Q , and superconductors [i[ . Vedral (H studied 
^ , the entanglement in many body systems at zero temperature using the second quantization formalism. Following his 

■ works, Oh and Kim Qlstudied the entanglement of two electron spins in a free electron gas, superconductivity Q 

| and the Kondo model [5[ at finite temperature using thermal Green's function methods. Recently, in the context of 

■ AdS/CFT duality @, entanglement entropy @, H, HI is used in studying the black hole entropy @, It was also 
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suggested (ill 1 1 21 ] that entanglement is an origin of mysterious dark energy. However, since entanglement entropy is 
a good measure of entanglement only for pure states, we need a reliable method to calculate entanglement for more 
general states. Thus, studying entanglement of quantum fields vacuum or thermal ground states is of fundamental 
importance (Tij in various physical fields. In (l5j . the temperature Green's function approach was used to investigate 
the quantum entanglement of two non-interacting spin 1 boson particles in thermal equilibrium. It is also well known 
that, in algebraic quantum field theory, due to the Reeh-Schlieder theorem (Trjj the vacuum for quantum fields violate 
Bell inequality and has quantum nonlocality [It], Oil ■ In Ref- [13; Gil entanglements of Bose and Fermi quasi- free 
(gaussian) states are extensively studied^while spatial entanglement of free thermal bosonic fields was studied by 
post-selecting certain momenta in Ref. |2C|. On the other hand, for decades the nature of the vacuum, concerning 
the Casimir energy 2l], 22], in a bounded space has been extensively studied using the Hadamard Green's function 



method [23(. In Ref. [24j we showed that this Hadamard Green's function is also useful to study entanglement of 
vacuum and used the method to investigate entanglement of massless scalar field vacuum with a Dirichlet bound- 
ary. Entanglement of the scalar field could be experimentally tested by the scheme with trapped ions [TJ, [36| or 



Bose-Einstein condensates [37[. The entanglement measures such as the purity [25| or the negativity has been also 
suggested [26] for two-mode Gaussian states. Since the scalar fields are continuous variables we need a separability 
criterion for continuous variables such as Gaussian states. To utilize positive partial transpose (PPT) criterion for 
Gaussian mode, we averaged the fields over two tiny boxes to make the vacuum state an effective two-mode Gaussian 
state. 

In this letter, we first generalize the methods proposed in Ref. [24j to arbitrary space-time dimension and use this 
Green's function approach to investigate the quantum entanglement of the vacuum for a massive thermal scalar field 
in the 1+1 dimensional cylindrical space-time. Studying the nature of quantum fields in a compact space is especially 
important in the string physics. Our methods could be used to study entanglement of continuum limit of harmonic 
chains in various situations. 

We start by introducing massive real Klein-Gordon fields, which are described by a Lagrangian density, 

A*.*) = ^(fla^V + mV), (1) 
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in (D + 1) dimensional space-time. The equal time commutation relations of the scalar field are 

[<P(x,t),<j>(x',t)} = 0,[7r(f,t) J 7r(f ,t)] = 0, 
[</>(af, t), 7r(a; , £)] = iS(x-x'), 



(2) 



where 7r(x, t) = dt4>(x,t) is a momentum operator for <f>(x,t). Consider a vector of the field and its momentum 
operator at two points x and x! at a given time t, i. e., £ = (</>(:?, i), 7r(x, t), 4>{x' , t), tt(x' , t)). The variance matrix for 
the vacuum or the thermal ground state of the scalar field (|0)) is defined as 



^ = -(0|{AC Q ,A^}|0), 



(3) 



where {A, B} = AB + BA and A£ Q = £ a - (0|£ Q |0) with (0|£ a |0) = in this paper. Following Ref. [2J], one can 
obtain the variance matrix 



V = 



a c 0' 

b d 

c a' 

d V 



(4) 



where a = (0\{(f>(x, t), <p(x , i)}|0}/2, 6 = (0|{7r(f,i),7r(f,i)}|0)/2, a' = (0\{<f>(x',t),<p(x',t)}\0}/2, and so on. Usually a 
and a' terms have a divergence which can be removed by subtracting a free-space Green's function Go from a Green's 
function G called the Hadamard's elementary function [27j ■ Thus we can define a regularized Green's function; 
Gr = G — Gq which is not divergent. 

In Ref. [2J] we showed that, once we know the Green's function we can calculate the components of the variance 
matrix V from it, i.e., 



c(x,x) = -{0\{<f>(x, t),(f)(x,t)}\0) = lim -G R {S, t;x,t'), 
d(x,x") = ^{0\{d t cb{x,t),d t (b(x',t)}\0) = lim / d t d t ,~G R (x,t ] x',t'). 



(5) 



Then, a — lim^/^j c(x, x'), b = lim^/^j d(x, x') and so on. 





FIG. 1: Cylindrical 1+1 dimensional space time with spatial circumstance R. we average fields over two tiny boxes B and B' 
centered at x and x' , respectively. 

Since 'if and only if ' separability test for infinite-mode states are unknown, we need to reduce the infinite-mode 
states to effective two-mode Gaussian states by defining collective operators. To define collective operators, we follow 
the approach in Ref. [T3 |. i.e., at a given time t we spatially average the field operator over two tiny spatial boxes 
B and B' centered at Xq and x'q, respectively (See Fig. 1). Defining collective operators in this manner is physically 
reasonable, because in a real situation probes always have finite spatial resolution. We do not average but integrate 
the momentum operators within the box, since momentum is additive. ( However, averaging the momentum gives 
the same final results in this work.) Thus the collective operators are 



(G) 
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Here J B d D yf(x + y) denotes an integration of f(x) over a box B centered at x with volume L D . Then, the 
commutation relations in Eq. (2j reduce to the canonical commutation relations for the collective operators [13] ] 



where 



n = 





Ep] = i 


flap, 


J 0" 


>J = 


' 1" 


J 


-1 



(7) 
(8) 



Similarly we define a variance matrix for the collective operators V a p = i(0|{S Q , H^O). 

Now we discuss separability of the vacuum. Using the fundamental theorem of calculus [28[, one can calculate the 
variance matrix V for the two- mode Gaussian states with those of V in the limit L — > (but still non-zero). For 
example, 



V 13 = |im i(0|{$(f,t),$(^',i)}|0) 



(9) 



1 



lim 

l^o 2L 2D 

1 



d u y / d^gm^x + y, t), <f>(x> + g, t)}\0) 



-(0\{4>(x,t),4>(x',t)}\0} = V 13 = c. 

This is always possible when the integrand of the second line is continuous. Similarly, 

1 



Bo L2D 2 L ™ 



d u y / d u g{0\{n(x + y,t)^(x' + y 1 ,t)}\0) = L 21J V, 



L 2D d 



(10) 
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(11) 



For L infinitesimally small the above equation gives an exact result, while for L -C 1 it provides us a good approx- 
imation. The separability criterion we use in this paper is PPT criterion [2^, [3(| for two-mode Gaussian states [3l| 
which is equivalent to 



F = S - (i +4detV) < 0, 

where E = detA + deti? — 2detG. By inserting the components of K into Eq. (|T2|) we obtain 

1 



L 2i, [a6 + a'b' - 2cd] - 4L 4U [aa'bb' - aa'd 



c 2 bb' 



c 2 d 2 ] 



(12) 



(13) 



Now the quantum separability can be tested for an arbitrary space-time using this quantity once the Hadamard's 
elementary function is known. Furthermore, if a = a' and b — b' one can immediately obtain an entanglement measure 
called the negativity using Eq. (9) in Ref. (32J; 



E(p) = max < 0, 



^(o-|c|)(6-|d|) 



1 



(14) 



Thus, we can easily determine how much a given state has entanglement, once we know the Green's function at 
least for the collective operators. This is very useful results since the Hadamard's function is well known for many 
situations. 

As an application, we test the quantum separability of the massive thermal scalar fields vacuum in 1+1 dimensional 
cylindrical space-time using the methods described above. This model is related to one-dimensional harmonic chain 
systems with periodic boundary condition [13l.l33l.l34l]. First let us choose two points x and x 1 (see Fig. 1). Using the 
method of images [35] one can obtain the thermal Hadamard Green's function for this field 



G 



-1 / 

- — In ( 16 [sin(7r(Au + i(3m) / R)sin(n(Av + ij3m)/ R)\ 

47T V 



(15) 
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where u = —x + t, v = x + 1 are the light-cone coordinates, and Am = — (x 1 — x) + (t' —t) = — Ax + At. This Green's 
function diverges as R — > oo; 

\ 2 / Aj. i A , •_ 

G . (16) 



lim G ~ - . 



/ln(i)^ ln(16(7rAt-7rAa; + i7rm (3) 2 (nAt + ttAx + in m (3) 2 ) _ 



4ir 



Interestingly, the Green's function is a function of m(3 = M. As mentioned above by subtracting Go from G one can 
define a regularized Green's function Gr = G — Gq, from which one can obtain 



(a, b, c, d) 
( 



2tt 



Trm j3csch(Z%£) \ -tt 1 



7? 



, —In 



16 tt 4 (Ax 2 + m 2 /3 2 



" 6i? 2 ' 8tt 



' 2tt 1 



n 2 CSC 2 (^) 



y v / y4i? 4 (cosC 2 -^) - cosh(^^)y 

using Eq. ([5]) after some tedious calculation. By inserting these components into Eq. (fT!?]) we obtain 

1 



i? 2 



F = — - + i 2 / 2 (Aa:, M) + i 4 / 4 (Ax, M), 



where 

/ 2 (Ax,M) = - 
and 

h(Ax,M) 



- 1 + 2 7T 2 Ax 2 + C0S(2 7T Ax)) CSC 2 (7T Ax) In ( 



(cos(2 7rA3:)-cosh(2 Af tt))' 
4tt 4 (Af 2 +Ai 2 ) 2 



ln(M7rcsc/i(M7r)) 



16 tt 2 Ax 2 
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(18) 



(19) 



576 tt 



Ax 2 — 7T 2 CSC(7T Ax) ) — 7T 



hi 



(cos(2 tt Ax) - cosh(2 M tt))" 
4 7r 4 (Af 2 + Ax 2 ) 2 



Here we have set R = 1, hence from now on all length scales are in units of R. 



16[ln(M tt csch(M tt))] 5 
(20) 




0.5 



FIG. 2: F as a function of Ax and M for L=0.01. 

Our numerical study shows that /i has a maximum /i(l/2,0) = 0.134 and / 2 has a maximum / 2 (l/2, 1.107) = 
0.0164. Hence, one can see that 

F < — - + 0.134 L 2 + 0.0164L 4 < (21) 
4 

for < L < 1/2. Therefore, the effective two-mode scalar field vacuum with the periodic boundary is PPT and, 
hence, separable, for all M and Ax. Fig. 2 shows F as a function of Ax and M for L = 0.01. 

Our results implies that changing mass, temperature, and periodic boundary condition does not induce additional 
entanglement to the possible entanglement of the zero temperature free space scalar fields vacuum. Our results 
also indicates that the entanglement of generic vacuum resides in the free space Green's function Go which usually 
gives rise to infinite quantities. Our approach provides a new and generic method to investigate entanglement of the 
quantum fields vacuum or thermal ground states in a various situation using the Green's function method. Since the 
Hadamard's functions are known for many cases our approach can be useful. 
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